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ABSTRACT 


A dynamic inventory control model, formulated in terms of 
discrete variable servomechanism theory, is examined. The 
behavior of the system is first analyzed by digital computer 
simulation. Both random and deterministic inputs are considered. 
Two double echelon models are then constructed and examined 
analytically utilizing transform theory. Finally, the proper- 
ties of forecasters of demand are analyzed for the cases where 
mean demand is constant and linear. Exponential smoothing 


techniques are compared with standard estimation procedures. 
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l. INTRODUCTION 

This paper will take for its starting point a dynamic 
inventory control model formulated in servomechanism theory 
by Vassian [6] and extended by Reilly [5] and Zehna [7]. The 
techniques of linear system analysis as developed by Howerd in 
[2] and the exponential smoothing/forecasting techniques of 
Brown [1] are heavily relied upon. The model as it is described 
in [7] will henceforth be referred to as the "basic model." 
Briefly, in the basic model we have demand as the input and 
inventory level, measured relative to some safety level, as the 
output. Inventory is replenished using a reorder rule which is 
based on a linear combination of past demand and inventory level 
feedback. Units of inventory are assumed to be issued instan- 
taneously and a fixed lead time is assumed for receipt of reple- 
nishment orders. In discrete variable servo theory notation the 
two equations which describe the basic model are 


1-1 = + 0 . d 
gS ants ees 


t t 
BI 
j™o j=o 


where I # on hand inventory relative to some safety level 
at time t. 
8 = reorder quantity computed at time t. 
t 


a set of constants which weight past demands. 


la } 


iB = a set of constants which weight past inventory 


levels. 


Ry 


Let D(z) =) 4, : , be the g-transform of the sequence’ {d, }. 


k=0 


Then, the above equations become, in transform notation, 


z ST+1 
(1-3) I(z) = %1(z) + & 8(z) = X(z) ~—s_ and 


de : 
4 


(1-4) 8(z) = A(z)X(z) + B(2)I(z). 


It is easily shown that 


rt 
aon 
“i - % = B(z) 


~~ 


Let 
_ grt 
ve A(z) «1 
(1-6) S(2) °C ————————— 
® gltl 


1-% “F B(z2) 

be called the iinctel function. Then.(1l-5) may be written 

I(z) = S(z)X(z). 
Similarly, 
(1-7) 0(z) = [A(z) + B(z) 8(z)] x(z) 

= $'(z) X(z) 

where §'(z) § [A(z) + B(z) S(2)] 
the basic model is summarized in flow graph form in Figure l. 


As Vassian suggests [6], we will let 


(1-8) Biz) = 20272) 
T+1 
Lez 


TACOW OLSVA ANT JO NVAOVId Hdvad Mold 


Lr FEN 





z) X 





sm 


It follows then that S(z) will have no finite poles. In that 
case, the denominator of S$(z) will be finite for all finite values 
of gz. 

The following additional relationships, developed in [7], 


are listed for ready reference: 


(1-9) E(z) = EL1(z)] = S(z) m(z) 


where 
m(z) = ELX(z) J 
T+1 
(gt 1 
(1-10) a¥(z) =z A(z) GZ) xy 
l-z 
T 
1-11) Ax = +) Xx 
( dA teLyr Tey tt+l+T-j . We require that 
jro 
t+1+T 
- ee. 
(1-12) RLA* J -) m(j)--0 as t-7@ 
tr LYr 7 
jzttl 


When expected demand does not vary with time (eLX, |] = a) we have 


the following additional relationships. 


using exponential smoothing. 


-k 
k=o 
; Y Xz 
(1-14) a(z) #7 X@. 
1 - Bae 
Get) Xe oT (T-1) a , a forecast of demand during lead time 


plus one unit of time. 
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Tt+l 
Y (T+1) zg X(z) 


1-662 


(1-16) A*(z) = 


a 
(1-17) A(z) a CS all 


T--1 , 
(1-2 )(1- 6 z) 





T+1 
T+ 
a (T+1) 3 ae) 
(1-18) s(z2) = ———-*—_ - + #§ — 
l- pa l-z. 
bh . 
De x T+1) 
COR, XS) 64 Bee 
(1-19) I(z) = 
l- ba, 1 «il 


a (T+) (1-2), | 


(1-20) 0(z) = $/(z) X(z) = X(z) 
; l - Zz. 
Invegting the transform in (1-20), we obtain 
t t-l 
(1-21) 0 (T+L) "Bk x @ (T+L) » 3 X +X 
= = 0 = 
tc t-k tel-k t 
k=o k=o 
and 
o = - 3 4 . 
(1-22) @ = @ (T+) Ix, - a +X, 
Inverting the transform in (1-19) yields 
t-T-1 T 
7 k 7 
(1-23) 1 =a (re) )’ BX -) Xe 
t 
t-T-l-k 
k=o k=o0 


ll 


Computing the expected value of I. we obtain 


teT-1 + 
Bie 
(1-24) E = a (T+1) ” eX m(t) - ) m(t) 
k=0 k=o0 
t-T-1 


k 
= a@ (T+1) a)" B - (T+l)a 


k=o 


= a (Tl) a= _- (Tl) a or, 


=H 
(1-25) E. = ea (T+1l) g° 


Thus, we have described the basic model and displayed the 
resulting mathematical relationships. In addition, for the 
case of constant mean demand we have derived explicit formulae 
for the expected on hand inventory and reorder quantities at 


time t when exponential smoothing is used to estimate the mean. 
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2. DIGITAL COMPUTER SIMULATION 


The simulation presented here was programmed in FORTRAN 
60 and run on the CDC 1604 serial number 1. The programs used 
are included in Appendix II in the order of discussion below. 

First, a deterministic simulation of the basic model with 


constant mean was made based on the relation 


t-T-1 T 
E = a(T+1) ) B* me) - ) met) 
k=o k=o0 


letting @ = .2, 8 = .8, andT = 4 In the first run m(t)=16 
and the results of 995 iterations are shown in Appendix I-A. 

It is seen, that to six decimal place accuracy, the expected 

inventory level reaches zero after 83 time periods. This is 

merely a simple verification of the analytical results shown 

in [7]. 

In the second deterministic run, we let m(t)=t in the basic 
relation. Since Reilly [5] has shown that first degree aso 
nential smoothing yields biased results when the demand is 
linear we would expect this bias to be reflected in the expected 


inventory. In the case at hand we have 


t-T-1 T 
E = a(T+l) » aX m(t-T-l-k) - ) m(t-k) 
k=o k=o 


Letting ® = .2, B = .8, m(t) = t andT = 4, 


t-5 4 
k 
E. =) (8) (t-5-k} - ) (tk) 
k=o0 


k=o 
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C—) t-5 t-5 k - & 
eo eee ke 
Bo=t) (.8) 5) (.8) - \"k(.8) -\'t +5" 
k=o ko k=o0 k=o0 k=0 
| 
t-4 | t- 
L- 8 } l- me) 
= t -5 Sy 


02 sl te 02 


t-4 t-5 
= 1-.8 = <2 t-4 8 
+ .8 sccmieany 022 Cae. eee ) aaa | -st + 10 


04 J 


= seuciaumt rm 45(.8)° _ RYE Oe 16(28) -35. 


Hence, lim E # <35 


t—> © 


We can see from the results shown in Appendix I-B that in 
fact E. = -35 after 78 iterations again to six decimal place 
accuracy. So we have verified the analytical conclusions in 
(5) and [7], i.e., that the model is asymptotically unbiased 
(no steady state error in the expected inventory level) when- 
ever asymptotically unbiased forecasters are used <cnd asympto- 
tically biased when correspondingly biased forecasters are used. 
However, as a periodic review inventory model, the basic model 
may require an impractical number of periods to converge to an 
acceptable level. Conceding that convergence can be accsierct 
by adjusting @ and T, we must face the fact that asymptotic 
unbiasedness may not be a practical measure of effectiveness for 


an inventory model. 
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The second simulation to be considered is stochastic in 
nature. Assume that demand is Poisson, i.e., 
sane ne . 
PLX =k] = —————-_ where 1 #)AD. 
t k! t 


For ease Of simulation a normal random number generator was 


empoyed taking advantage of the fact that for large A [3] 


k: 





The 1000 random numbers for the first two runs are listed in 
Appendix III-A, while those used in runs three and four are 
shown in Appendices III-B and III-C respectively. 


In the first run A = 16 in the basic relation 


t-r-1 T 
ye Dx 

ee Keet-1-k “L"t-k * 
k=o k=o0 


ae thru ee are listed from right to left in Appendix I-C. 


It is seen that there is some tendency toward stability through 
Las but after that time period we observe the effects of the 
random demand on the inventory level. It is now clear that 
carrying expected values through the model does not adequately 
describe the true situation. 

In the second stochastic,run, all parameters and variables 


were unchanged from the first run except that T was increased 


to T = 10. We see the resulting inventory values in Appendix I-D. 
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As expected, some increase in the variation of r. is observed, 
i.e., the variance of the random variable I. is a function of 
the lead time of tie basic inodel. However, we also observe 
that Ee is relatively insensitive to T in the sense that the 
order of magnitude of the change in T is not carried over to 
tT to an appreciable degree. 

In the third stochastic run, all parameters and variables 
were the same as in the first run except that de 100. Natur- 
ally, to preserve our Normal approximation, the Poisson distri- 
bution variance was also increased to 100. Note that the random 
numbers generated and listed in Appendix IlI-Dreflect this 
change. Now we observe tiie resulting values of I shown in 
Appendix I-E. Analogously we see that the variance of _ 
depends on the variance of - but is relatively insensitive. 
The order of magnitude of the increase in A is not reflected 
in Is 

In the final run Lv was allowed to vary with time, speci- 
fically d. = 16+ +t. From the analytical development we are 
led to expect bias since our forecast is biased for non-con- 
stant mean demand. In Appendix III-C observe the new set of 
random numbers generated with mean A = 16 + t and variance 
Ais = 16 +t. The expected bias is apparent from an inspection 
of the Appendix I-F. Most values of I, are negative. Again 
observe the relative insensitivity of - to increasing the 


variance of the demand. 
In summary, the computer simulations have confirmed the 


theoretical results as expected. , However, in doing so, some 
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limitations of the basic model have been suggested. Asymptotic 
unbiasedness as a criterion appears to be week and probably 
impractical even when expected values of the random variables 

SS and I, are considered. When {x, } is a sequence of realizations 
of a random variable from a probability distribution, the sequence 
{z,} certainly does not appear to approach zero, Vassian's [6] 
apparent claim to the contrary. Zehna [7] clarified the theory 

in this area but the results of the simulation dramatically 
emphasize the inadequacy of the basic model when only expectations 
are considered. The choice of a safety level, or other wandeeeial 
decisions, cannot be made on the basis of eli, agdone but var(1,] 


must be accounted for also. 


Lf 


3. DOUBLE ECHELON MODELS 


It is pOssibie to formulate multiple echelon models in 
servo theory along the lines suggested by Howard [2]. It is 
desired to test the basic model of the previous sections in 
echelon. First, we construct a double echelon model using 
three basic models, two at the retail level and one at the 
wholesale supporting stock point position. Assume that con- 
sumer demand with constant mean values a, and = is the system 
input at the retail level. In addition to the inventory level 
outputs at the retail stock points we have the reorder quentities 
8 (t) and 9,6 serving as demand inputs to the wholesale stock 
point. There is a lag time which is constant for each retail 
stock point of T, 


by the wholesale stock point and receipt by the retail points. 


and 5. respectively between the time of issue 


Assume that the wholesale stock point has a sufficiently high 
safety level so that it is always able to fill orders from the 
retail activities. There are to be no interactions between the 
retail stock points. All three stock points are to use a single 
exponential smoothing forecaster which is of course asymptotically 
unbiased for constant mean demand. The lead time, tT) for the 
wholesale stock point is the delivery time from the manufacturer. 
Assume that all three stock points use the same value of the 
smoothing constant @. Hence this double echelon model describes, 
in a simple way, a common situation involving two echclons of 
inventory. The use of a common inventory control doctrine and 


the assumption of no interactions at the retail level are considered 
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to be reasonable in the light of the author's experience. It 
is then of interest to determine the possible effects of such 
an overall doctrine in terms of this model. 

In the framework of the multiple input-output system in 
[4] we have a system with two inputs X (2) ; oe and three 


outputs I, (2), I, (z) and 1,(z)- In matrix notation this becomes 


f x 
h (2) 0 h (z) 
13 
[1 (2), 1,(2), 1, (2) ] = [x, (2), x, (z)] | 
T T | 
0 h (z) h (2) | 
i. 2 ae 
where, for the present case, 
ac 
T, +1 i 
T Zz (1; Hi) Il-z 
Rs ES) na rs 
ll . l-62z lez 
T, +1 
+1 
5 z* oT, +1) 4_,72 
and h (z) =S (2) 2 ———-_ —--- —*—_.. 
22 “ 1-Bz ie 
= nd 
Hence, T(z) S, (z) X, (z) a 


aCe) ~ S,{) X,(z) : 


This double echelon model is graphed in Figure 2. Nov 


I (z) = S,(z) [@_(z) + 98 (z)] 
3 4 1 2 
since the input at the wholesale stocking point is 


8, (z) + de 
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TA0OH NOTHHOd ATENOG 
ALVSTIOHM/TIVIdd 20 WVdOVIG Hdvad MoTtd 


co AYUNIA 





ITS ATOM UES 





From (1-7) we have 


De > S ,() (La, (z) + B, (z) S(z)] X, (z) + La, (2) +B, (2)S,(2)] 


x, (2) 3 


Hence, 


T 
i = S ,(z) La, (z) r B (z) S,(z)] and 


T 
—" = — LA, (2) + B (2) S,(z) ] : 


Recall that we ate asswalag 


m(t) = ELX, (t) J =a 


1 and m,(t) = eLXx, (t) ] = fe for t 20. 


Hence, 


a 


m, (2) oe and m, (z) & = “ 


Consider now 


t k& t-1, 
8 (t) =o (4) \'8 — - @ (T+) > BX yy? XO) 


k=o k=o 
so that a Awan 
9 (2) = eee + * X.(z) 
for i=l, 2. 
t t-1 
Ele (t)] = a(T +1) a, i - &(T +1) _ \" B+ “a 
k=o k=o 
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(B'=1) (B=) 
“ ah Sg S ain eae +a 


BLe (t) J o(T +1) a. 


B-1 


t t-l 
a [1-(r +1) (8-1) + (T,+1) (8 -1)) 


t tol 
a [l-(7 +1) B + (T +1) 8B j 
1h 1 i 


-1 
a Cie(r +1) (BR - 8°)] 
i 


1 
(1 - B)] 


t 
a (1+(T +1) B 
4 1 


t-l 
a [1+(T +1) a B i = 
i 1 


This demonstrates that the expected demand at the wholesale 
stock point, which is the sum of the order quantities from the 


retail stock points [8 (t) + 9, (t) J], is not a constant despite 


the constant mean demand inputs, m, (t) and ae to the system 
as a whole. It is clear then that the inventory control doctrine 
at one echelon does affect the stability of the system. In this 
case, a cOnstant mean demand is transmitted as an exponential 
function because of the exponential smoothing forecasting system 
emp loyed. 

Now it is of interest tO examine the behavior of the whole- 


sale inventory level under these conditions. Consider 


22 


a La, (2) Bs B, (z) S (2) ] X, (2) 


which is the first component of 1, (2) shown above. By assumption 


_ 


T.+1 T +1 
I aie 
3 l1-B2 lez 


[a (2) + B (2) S,(2)] X, (2) 











Pe 
@(T,+1) (1-z) _ 8 a(T +1) (1-2) ae 
T +1 T + 
(1-z + )(1-Bz) (lez 1) (1-8 2z) 
a 
ts 
l-z 
a(T +1) (l-z) + (1-8 z) a 
Cle a LO 
1-68 z l-z 
Hence, 
5 (ala, (2) + By(z) 8,(2)] X (2) = S,(2) S) (2) X, (2) 
4 Tt T,+1 
rst 2 = 3 a(T +1) (1-2) +(1-Bz) 
: ee ake 1-Bz 
a 
— 
l-z 
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3 
a(T,+1)z ? Z 
in | eee eC from [7, p. 15], 
z—>il 1-68 z 1-2 


and obviously 


@(Dy+1) (lazy (1-8 oz) 


Lim ie <4 
z—>1 1-6 z 
Therefore, 
lim (l-z) S (z) = 0 
2-1 3 
and since 


B (2) =8,.(2) [s. (2) X(z) + S$ (z) X(z)] , 


it follows that 


lim (l-z) E (z) = 0 ‘ 
z2->l1 3 
Based on a well known theoren leds, », 210], the result above 


implies that E79 as t-°, 


At first this may seem to be unexpected until one realizes 
that asymptotic unbiasedness is not a necessary condition for 
the steady state error to vanish. Recall that, in transfer 


notation, 


(lezJE (2) = S (2)8{() a + 828) a2, 
#8, a, + 8,3) 2, 


Thus with constant mean demand, the most that is required for 


asymptotic unbiasedness is that lim a2) = 0 and that 
zZ--l 


Lim S53 (zs) be a constant. Since S (2) depends only on A(z) 
Zz a) 
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and B(z), it is relatively easy to choose A(z) and B(z) to yield 
the limit zero. S, (2) here is a linear combination of the 
asymptotic unbiased retail level transfer function and hence will 
be a constant in the limit under these conditions. Therefore, 
the results above are indeed natural under the circumstances. If 
initial expected demand is not a constant (m(z) # a /a-z) then 


ie and (2) for some constant c may not be sufficient 


conditions for no steady state error. Again this emphasizes the 
weakness Of asymptotic unbiasedness as a criterion. 

Inventory managers have a keen interest in the number of 
units in the pipeline (units issued from one stock point but 
not reflected on the records of another). Let us see what this 
double echelon model reveals about the problem. Remember it was 
assumed that issues are made instantaneously by the wholesale 
stock point but sae andr ot? time periods elapse before the 
inventory is recorded at the respective retail points. We are 
interested in the expected number of units in the double echelon 
system which are not shown on the records of any stock point as 
inventory. 

The total number of units on order at time t by a retail 


stock point can be written 


t=] 
a(T +1 

Q SQ," ) > meatier 

jal j=l k=0 

Hence, 
it t-j 
k 
BLQ J a (t+) ue 
j=l k=o 


25 


T t-j 








eel 


T 
eCQ) = a(t+l) a\’ 7 ok sare) a Y (B=) 
Le is 
j=l  k=o j=l B- 1 
T ; T ~i4] 
t-j+l t-j+ 
= -a(rti) \(B-1) = at(ttl) -a(r#1) ) 8 
j=l j=l 
gay . t-T+2 
= aT (T+L) -a(T+)) ca) + a(T+1) 20 
B-1l Be-l 
(T-+1) cel t-TH2. 
= aT(T+l) + a B -f | 
t-T+2 
a(T+1) 8B T-l 
= aT(T+1) + (B =1) 
a 
t-T+2 T-l 
a(T+1) 8 (1-8) 
= aT(T+l) - : 
oO 
Then, lim EQ ] = aT(T+)) 
t—>© : 
for any retail point with ELX, J = 4 


and so 


Lim ELQ J 
t -@ : 


aT (T +1) + aT (T +1 
TT, +O tat a+ D 


for the present double echelon model. 
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This is a natural result and the steady state condition is 
a direct consequence Once again of the asynptotically unbiased 
forecasting. Note that the expected amount in the pipeline is 
a function of the square ofvleal time wn.ch re-emphasizes the 
possible savings resulting from a reduction in the lead time in 
any veal situation. Then it appears that models of this type 
are useful for further study of pipeline inventory. 

An inventory control point (ICP) is an office whose function 
is to manage a so-called commodity area which usually consists of 
a large number of stock items. A substantial portion of its 
effort involves the maintenance of optimal inventory levels at 
a system Of stock points. The stock points report changes in 
their inventory levels and forward replenishment orders to the 
ICP periodically. The ICP takes centralized procurement action. 
One point of great interest is the correlation between inventory 
at the ICP and the inventory as it really exists at the stock 
points. It is possible for incompatible inventory control doc- 
trines at the various stock points as well as different lead 
times, to generate a record of demand and inventory at the ICP 
which differs greatly from the actual demand input to the system 
and the inventory at the stock points. <A very simple double 
echelon model which operates on the basis given above can be 
analyzed. 

Assume there are two stock points operated as basic models 


with constant mean demand. Again the mean demands ) and ay are 


not necessarily the same nor are the lead times Tt and T.: There 


oy, 


are no interactions between the stock points. The reorder 


quantities generated at the two stock points, 9 Ce) and © (t), 


are transmitted to an inventory control point wiose operation 
differs from that of the basic model. It is assumed that, 
unlike the stock points, the control point has no physical 
inventory. It merely keeps a record of system inventory carried 
at the stock points based on the inputs it receives. Besides 
the two order quantities, the ICP receives a report of the 
inventory level at each stock point at each time period. Material 
ordered by the ICP is delivered direct to the stock points so 
that the lead time T, is a part of T, and ae The control point 
does not feed back its recorded inventory level but utilizes the 
values of I(t) and T(t) to determine order quantities 2A) 
and update I,(t). The basic equations which describe the oper- 


ation of the ICP are therefore 


= er “1) € 9 (t-— - (0 ) 
1, (t) L(t 1) + z(t 1) st t T,) [ 6t) + s(t) J 


t t 
,(t) “Sa, [0 (t-3) + 8 (e-3 J +) B, [rT  (t-J) + 1,(t-j) J 
j=o j=0 


and in transform notation 

T,+1 
I (z) =z 8 (z) - 6 (z) - 8 (z) +z C1.(z) + I (2)) 
3 3 iL 2 1 2 


8 (z) =A (z) [6 (z) + © (z))] +B (z) [I (z) + I (2)] : 
3 3 1 2 3 1 2 
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The flow graph for this double echelon model is shown in 
Figure 3. The inventory contro! » «at's system is graphed in 
Figure 4. 

Now let us examine the behavior of the ICP inventory record 


in relation to the actual inventories. 
T3tl 
z [I_(z}+I (z)] +z A (z)L[®_(z)+0 (z) J 
1 2 3 1 2 


T(z) 


+ B(z)[I (z)+I 1. [8 (z)+9 (z)] 
1 2 1 2 


TQ.tl 
3 
alS, (z)X, (2) + S$ (2)x, (2) J +z {aot jen core; (0 


+ B(z)Ls_(2)X, (2) + $5 (20%, (21) 084 (a)x (248° DK, (2) 


T,¢1 
3 
[as (2) +2 A, (2) (2) + B (2)8, (2)-S (2) x, (2) 


T 41 
+[zs (z) +z A (z)S‘(z) + B(z)S (z)-S‘(z)] X (z) 
2 3 2 2 2 2 


S (z}X (z) +S (z)X (z) 
13 1 23 2 


where 
eal 
‘ 3. t t 
S (z) = 28.(z) +2 & (€2)S° (z) + B(z)S (z) - S (2) and 
13 t 30° 1 1 1 
+1 
a2 


S,4(2) = z8, (2) +z AtOe oe + B(2)S_ (2) - (2) 
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INIOd ‘TOUINOD AYOINAANT FHL 10 HVYOVId Havad Holts 


& AYNIDA 





LNICd TFOULNOD ALOLNFANI 


"69 + E)'O 


Ya. | 


Since 


Lim S (2) = 0, lim A (2) = 1, 


zZz—l Z—> 


1 


T +1 , 
3 





lim s‘(z) =1, lim B(z) =- 
z—->1 } zl 


it follows that 


lim S (z) = 1-120. 
zg—>1 13 


Similarly, 


lin § (2) =0. 
z—->1 23 


But, 


E (z) =S (z}) m (z) +S (2) m (2) 
3 13 1 23 2 


a a 
1 2 

=S i)" +S 7) 
13 l-z 23 l-z 





and 


(l-z) E = a, S14 (2) + . 2G 


Therefore, 


lim (l-z} E (z) = 0, 
Zo 3 


which implies that 


lim £E (t) = 0. 
t-> © 3 


Again there is no steady state error. Basically the results here 
parallel those in the first double echelon model. In effect, 


major difficulties can exist in the type of systan discussed 


re 


above but the "no steady state error" criterion does not really 
measure the behavior of the output function adequate:y. It 
appears feasible, however, to formulate multiple echelon servo- 
mechanism inventory models which include ICP type operations in 
Order to gain further insight into the Operation of such a 


complex secon. 
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4. DEMAND FORECASTING 

Reilly [5] and Zehna [7] apply exponential smoothing to 
define forecastors in the cases where mean demand is constant 
and also when it is linear. Moreover, in terms of the basic 
model presented here, Zehna [7] has demonstrated that the 
expected inventory level converges to zero whenever any asymp- 
totically unbiased forecastor is used. For example, in the 
constant and linear models, if least squares estimators are 
used to estimate the parameters and these are substituted in 
the forecastor, then, being unbiased, the forecastor is, a 
fortiori, asymptotically unbiased and the inventory level will 
converge to zero. The justification, then for using exponential 
smoothing lies in the advantages which relate to digital com- 
puter storage and mathematical tractability in the servomechanism 
model. See Brown [1] for a complete treatment of the advantages 
of exponential smoothing. However, these advantages are external 
to the present requirement of unbiasedness in the limit. 

Under the model m(t) = bt, Zehna [7], following Reilly (5], 


[2] (ou) 
t 


(x) ] where S$ (X) represents the 


b ——=\(isecxyuens 
uses b. 5 mes 2 


sequence i doubly exponentially smoothed. It follows that b 
t 


is an asymptotically unbiased estimator for b sv sit thea fora 
caster t+T+1 
——? A . 
A* = y b J 


is an asymptotically unbiased estimator for total demand over 
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the next lead time. For this model, another possibility suggests 


itself quite naturally. Since, 


e[s (x)] = s (2) Soe - 
t t oO 


ELS, (x) - x(t)] = - bid 
cvs 


and 


am 
lim EL Z (X(t) - S,(x))] = b 


t—©® 


Hence, the estimator 


“A 


b s us 
2 7 (x(t) - 8,09) 


is also asymptotically unbiased for the parameter b. Using this 


to define a forecaster 


ttT+1 
AX = ‘- a | 
teT+1 ee it follows that 
jettl 
t+T+1 t+I+1 t+T+1 
ae 3 r G7 ae. 
elas Jee) %5)* Fo by =) mp 
t+T+1 
j=tt1 j=ttl j=t+l1 


so that A* 1 defines an asymptotically unbiased estimator for 
t+1+T 


demand. Again by Zehna [7], lim zE* O when this forecaster 
t-30 


is used. Since this basic result is the same as that obtained 
when double exponential smoothing is used, other means will have 
to be used for comparing the two. 


First consider, 
{2 [s (x) “S. eo { = 7" [s.(8) - S_ (8) 
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i 22 At 
ei ([s@es. Glib=s 212 cape Eat awwecnane 
3 = L 8 t a 


—_—. 


; bow ce t-1, 4 t(ttl)b gt } 
= § (1-8 + Roles (1-t)-tB J+ FC 1-6°l- f 


(4-1) . -< ee = 


t(t+1)b apt 


+1 
a {bB . bB t(t+1)b age | 
at ae . a 


2.t 
t t(ttl)b a Bs 
=<b-bpB - ————— P 
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On the other hand 


£ (Stas) - 5,00) [oe = 8,001 
B tc @ B&B t 


(4-2) “ - : t-1 
= —<bt - § (1-8 )+ = [1-B (1-t) - tB 1 
Bien c - 


€ t-1_ 
= ‘bt - bt + btB + Be [1-8 (1-t)-t8 1| 


WW 


, 


(5 ie t , 
=~ 2-4bth + bB ri-p + 5 Gar 4 
Bt i 


Be i) ee Gee 
B 04 eg 
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bts 


Mes [x(t) - s co) | ELI Ae ye (B-1)+ best J 
8 t B ‘og 


1 
: 2 bB ip) ] = —% BB - pp ] 
B w B a rs 


ct 
= b - bp o 


Obviously as t-&° (4-1) converges to b more rapidly than (4-2) 


does and so is preferable. 
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5. CONCLUSIONS AND RECOMMENDATIONS 


The main conclusions reached here are that expected 
inventory levels and asymptotic unbiasedness appear to be 
rather weak criteria for the models considered in this paper. 
It would be highly desirable to be able to relate costs and 
effectiveness to these models. The inventory variance should 
be the basis for any discussion of effectiveness (expected per- 
centage of demands issued) and the quality of the estimator. 
Average cost depends on inventory carrying costs, order costs 
and shortage costs mainly. It is important that the decision 
criteria be related to these costs. 

Many possible courses of action involving the ramifications 
of the basic model exist. It is recommended that the following 
two extensions be attempted first to determine whether or not 
this type of model has sufficient power to be useful: 

(1) Calculate the Varlt J for several types of random 
demand and forecasting schemes. Explore the methods of deter- 
mining a safety level to provide a given degree of protection 
against shortages. 

(2) Attempt to define a new criterion for acceptable 
inventory level behavior which can be related to the commonly 
considered costs discussed above. Some insight into the pro- 
blem might be gained by relating the lack of steady state error 
criterion to costs. Specifically, cost would seem to depend on 


the amount the inventory level varies from the norm, the relative 
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values of the cost components, the number of times I. changes 
direction and the rate of convergence of ae 

For this model to be useful it is necessary that the 
reorder policy be optimal in the sense of minimum cost. After 
a minimum cost criterion is applied to the model, the stochastic 


process i t = 0,1,2,+-+} must be intimately understood. 
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aS RS ED eG Te SL -53,081 -72,465 -64, 
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-49.795 -53,436 -75.,549 -41,439 Sloot 
-19.681  -3$7.3845 2,924 °+».-4,461  £#-13,969 
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~ PROGRAM MEAN 
DIMENSION &(1000),e:(1000) 


_ CO THIS PROGRAM SIMULATES MILEL WITH MEAN DEMAND INPUT 


DO 1 181,1000 
PR = | 
90 4 [56,1009 —_ — eee 
SUM=0.0 
Me T=“5 
Weve, J=i5 <<< 
TEMPLER( J -4-J)* (6. B)aw( dad) 
2 SUM=SUM+TEMF1L 
TEMPL=SUM 
SUM=0.0 
DO 3 ae, 
_» TEMP2ER(iri-y) 
5 SUMSESUM+TEMF? 
Temegd=eSum 
4 (1) 21 EMPL: 1 EM 
PRINT 5 
5 FORMAICLLHMNEAN VALVES) 
6 FORMAT (6K,29F195. 6) 
ENT 
END 
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PRUGRAN. PO.1sos 


= eee i ee ee ee a ee ee —- 


C THIS PROGRAM SIMULAIES BASIC MODEL WITH RANDOM DEMAND 


Te ananaaniemeinanesaaiiad 


ee eS 


—DIMENSTON 


1 C0 re Oiled OL ie 





se epg geen a tn mt 





COMMON TRANDOM-.RANDOM 


TRANDOM25574655 
[= i. 1h 

eat 
bx | JSP? 


5 il ala CE 


EP Fe 


Sec ARIS 


a FR —P  P ID  OEH ENG AF PEED SE 


oa ee re eer 


res te SR SE 


_ 8 FORMAT(6X.5F10 4) 


Do ¢ Is42,1000 
SUM=0.0 
hace 
DO a J=i1,™M 
TEMP1=2,2 | 
2 SUM=SUM+TEMF1 
TEMP1=SUM 
Si eae C 
pew Js 1 yee 
TEMP@=X(f-JUr1) 
S SUM=SUM+7TEMP2 
TEMP2=SUM 


eT COTS S TEMP ia ee 


PRINT 5 


3 _ FORMAT (C2SHRANIUM NUMBERS GENERAIED)_. 
RSCG) 


PRINTS xs. i 


me ew ee 


NORMAL (16.0,4.0+X(1))~ 


/— RN oe a se i eG a i 


15) 1 ee ee CL = L0H oe 


rR ete Ae mi 


AO A TD AY 


POL RD aa A A ee ee EE 8 en A otrett 


G6 FORMAI{( 6X; oF 20 oe 


PRINT 7 


7 FORMATCL7HINVEN TORY ON HAND) 


PRINT 8, (101), 1222,2000) 


END 
SUBROUTINE 


ICAN, 
TRANDOMs 


2 Ys RANDOM 


RANDOMSY/67108864.0 


RETURN 
END 
SUBROUTINE 
COMMON 
SUM=0.0 
pO 2 get ue 
CA La RANVAR 

L SUM=SUM*RANDOM 
a IM = 6 ae 
Cakes f 
REI URN 
EN i 


NOK I 


AR 
COMMON JFANDOM.RANDOM. Ou 
CON(K1=67108864) 
fRANMDOM* S125 

IF CTRANDOM-~67108868)2,24,1 ae 
1 ENO(G) LDAtLRANDOM) »OVI(KL),$TIC1RANDOM). 


ecu EO, S64. cE 


=== 


ALC A,8,0) — 
[RANDOM.» RANDOM 


62 


3 


aheieeniate ame oieneeneet — 


C TH 





5 FORMATC24HRANDUM NUMBERS GENERATED). 


COMMON JTRANDOM.RANDOM | 


aa a a hy 1 ae 


—_——— eee ee A i 
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PROGRAM Of tain 
IS PROGRAM SIMULATES BASIC MODEL WITH 2ANDOM DEMAND 
EWE SO Nes G2 U0 Cr OO et ny ele 
COMMON TRANDGM-RANDOM 
TRANDOMs 3574655 
Uae seat 1 UDG en. — 
CeO) eG 1 Oe Gee. OO (1.).) 
KC L)=X¢(1) 





Oe tere 6,1 0 00 
SUM=0.0 
le eames A: 
DOwicmes) =o. 
TEMP1=(,.8)¥*( J-4)¥X( 1-4-9) 


te 
ee, we ————— 


2 SUM=SUM+TEMP 4 


ee -< ~ earn ennteentenensimnetiaatiindinemnemnernetenee none tiie ether ee ademas 5 ~e 


Tene SuM 
SUM=9.9 
06 Siri. 
TEMP2=X(i-u+1) 
S SUM=SUM+TEMP2 
_TEMP2=SUM . AE eee eee ee 
4 T(1)=TEMPL-1EM. 2 
PRINT 5 





Pe hiwowec AC 1), i=, 000) 
Omir URMa i coAs oF 10. 3) 
PRN a 
7, FORMATCA7HINVENTORY ON HAND) 
enema a | Gl jgchk=6. 20009 
oO cw wok ot LS) 
END 
SUBROUTINE RKANVAR 





CON(K1=67108869 ) - 
TRANDOM=IRANDGW 23125 
TF CTRANDOM-67198865)2,2,14 


ee 





ce ENG(O)»sLDACIXANDOM) ,OVI(K1)sSTOCIRANDOM) 


2@ Y=TRANDOM 
RANDOM=Y/67198664.0 
RETURN | 
END 

Mole. Oe LINE NORMAL CA, Bo) 
COMMON JTRANDOM, RANDOM 
SUM=0.0 
Homes he. 
Cheer AN YAR 

1 SUM=SUM+RANDCM 

= Sik MissOyel) 

C24 tore 
RETURN 
END 
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PROGRAM POTSSO™ 


GeTHTS PROGRAM SIMU Ee. BASIC MODEL WITH RANDOM DEMAND 


—PIMENSTON.1¢2000),X(20000,1x(1000)_ 
COMMON TRANDOM-RANDO4 
TRANDOM=$574699 


7=1 
WwW  Oealree 
eae NOPMALT Il xe Ge 
LXUigent 1? 


VESOR TCH sree ‘ we ee ee 


—— = a ee 


Do iodo 1000 = =~ See re 


ae xX CL 


yO 4 2ouee00 
| SUM=0-0 
ee 
DO 2 ORR , 
TEMP12(.8) 8 (Ur) ¥x( 1-4-4) 


oe AE eh RD TET + 


2 SUMESUM+ TEMPS _ | 2. 
TEMPL=SUM 
SUM=0.0 | 
DO 5 J aie peg <_ a. = oe See SP —ad 


= ~TeMP2eX(I-utl) 
3 SUMSSUM+TEMPe? 
4 TC(L)STEMPL- 1 EMF 2 
a1 Nias 


PRINT Te ae is ee re Mose 
6 FORMAT (6X,5F 10:9) 
PRINT 7 - 
7 FORMAT(LZHINVENTORY UN rAND) 
PRINT as oii ap ra cre) 8) 


8 FORMAT(6X,2FI0.9) ee a 


END 
SUBROUTINE RANVAR 
CON(K12=67106864 ) 
TRANDOMSIRANDOM# S129 
Io CIRANDOM-67108869)2,4,4 . — 
ENO(O) LDACIRANDOM) -DVICKL), STACI RANDOM) 
VeTRANDOM 
—RANDOMEY/67108864%, 0 
RETURN 
END 
SUBROUTINE NORMAL (A, 8,0) 
COMMON TRANDOM- RANDOM 
SUM=0.0 
CALL RANVAR 
4 SUMsSUM+#RANDOM 
Xa SUM~6,0 
C2x*Bra 
RETURN 
END 
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